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Abstract 

The perturbative QCD corrections to the nonleptonic decay rate of the b quark are 
discussed. By considering the limit where the final charmed quarks are slow, it is argued 
that the coefficients of the a 2 s terms, corresponding to next-to-next-to-leading order in the 
standard renormalization group expansion in hi(mw/mb), are naturally large. The large 
coefficients arise from the final-state gluon exchange between quarks and are associated 
with the region of rather low momenta, which may further enhance the phenomenological 
significance of these terms. 



1 Introduction 



The problem of the nonleptonic decays of the B mesons is well known and is well documented 
in the literature (see e.g. in [|IJ, 0): the experimentally measured semileptonic branching ratio 
B(B — > / v X) is too low, and leaves unaccounted in the existing theoretical calculations an 
enhancement of about 20% or so of the nonleptonic decay rate of the B mesons relative to 
their semileptonic decays. One effect: a relative enhancement by possibly as much as 30% 
of the sub-dominant decay b — > ccs due to the mass of the charm quark in the 0(a s ) QCD 
correction (which effect is contained in the results of Ref. 0, and which has been rediscovered 
more recently® 0), is not entirely sufficient for the overall enhancement of the total 
nonleptonic rate and also does not look to be supported by the data. Indeed a solution of 
the problem of the semileptonic branching ratio of the B mesons by an enhancement of the 
decay b — > ccs would require an average yield of charmed quarks per B decay of about 1.3, 
whereas the latest experiments result from CLEO for this number i S n c = 1.15 ± 0.044, 
which does not support a relative enhancement of this decay mode with respect to the 
dominant one b — > cud. Thus the experimental situation suggests a further theoretical study 
of the nonleptonic b decays beyond the thus far considered terms. 

A theoretical study of the ratio of nonleptonic to semileptonic decay rates of the B 
mesons involves two ingredients: perturbative QCD corrections to decays of a b quark and 
the nonperturbative effects related to the fact that the b quark decays being confined in a 
hadron. The latter effects, when considered within the expansion in the inverse mass of the 
b quark, are estimated'ED to be quite small: of the order of few percent. Thus the required 
by the data enhancement of the nonleptonic decay is presumed to be associated with the 
perturbative QCD effects in the decay of the b quark. The standard way of analyzing the 
perturbative radiative corrections in the nonleptonic decays is through the renormalization 
group (RG) summation of the leading log terms and the first next-to-leading termsS in 
the parameter L = lii(mw / mb) ■ F° r the semileptonic decays the logarithmic dependence on 
mw/mb is absent in all orders due to the weak current conservation at momenta larger than 
rrib, thus the correction is calculated by the standard perturbative technique, and a complete 
expression in the first order in a s is available both for the total rate® and for the lepton 
spectrum 113 . In reality however the parameter L pa 2.8 is not large, and non-logarithmic 
terms may well compete with the logarithmic ones. This behavior is already seen from the 
known expression for the logarithmic terms: when expanded up to the order a 2 s the result 
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of Ref . H for the rate of decays with single final charmed quark takes the form 



r(6 — > cud) + F(b — > cus) a s a 



1 + — + 



3 T(b — > cev) 7r 



7T 
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where, in terms of notation of Ref.[|J, c(a) = 022(0) — 012(a). The behavior of the function 
c(a) is known explicitly® and is quite weak: c(0) = 19/2, c(l) = 6, and c(m 2 /m 2 ) pa 9.0 for 
the realistic mass ratio m c /mb pa 0.3. One can see that the term with the single logarithm 
L contributes about two thirds of that with L 2 in the term quadratic in a s . Under such 
circumstances the RG summation of the terms with powers of L does not look satisfactory 
for numerical estimates of the QCD effects, at least at the so far considered level of the first 
next-to-leading order terms. 

An alternative way of evaluating the QCD corrections would be to use the straight ex- 
pansion in the QCD coupling at least to the order a 2 , which amounts to a calculation of 
the non- logarithmic term in the ratio (0) in the order a 2 . In terms of the RG summation 
this corresponds to calculating the next-to-next-to-leading terms, which can be used for con- 
structing the corresponding full RG expression in that high order. Unfortunately an explicit 
calculation of the a 2 terms poses a quite complicated technical problem, which however, 
perhaps, can eventually be solved. The purpose of this paper is to present arguments that 
the non-logarithmic terms in the order a 2 may naturally be large and in fact may compete 
by their numerical significance with the known logarithmic terms. In other words the natural 
parameter for these terms is a 2 in the case of b — > cud(s) decays and even 7i 2 a 2 in the case 
of the decay b — > ccs rather than the (a s /7r) 2 that is usual for the radiative effects. To see 
this the limit of small velocity of the c quark in the b — > cud(s) decays and of small velocity 
of both charmed quarks in the decay b — > ccs is considered (slow charm limit). The realistic 
kinematics of the b quark decay is knownlilt to be not too far from this limit, however it 
is far enough to prevent from obtaining a reliable quantitative estimates in the discussed 
problem of the real b quark decaysQ The reason for the large coefficients, arising in the limit 
of slow charm is the Coulomb-like gluonic interaction of the quarks in the final state. The 
presence of this interaction makes the nonleptonic b quark decays fundamentally different 
from the semileptonic ones at least in the slow charm limit. The QCD corrections to the 
semileptonic decay are finite in this limit and are determined by exchange of virtual gluons 

1 It can be also noticed that for the decay b — > ccs the inverse of the velocity v of either of the charmed 
quarks in the center of mass of the cc system, averaged over the phase space in the decay, is ~ 2. Thus 

v^ 1 is not a much worse parameter than L. However, as will be seen in the real kinematics, an expansion in 
v of the a 2 s terms lacks convergence similarly to the expansion in L. 
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with the scale of momenta between m c and m& Hi On the contrary, the QCD corrections 
to the nonleptonic decay rate are singular in the limit of zero velocity of the final charmed 
quark already in the coefficients of the a 2 terms. 

The set of essential graphs, which need to be calculated for a complete 0(a 2 ) QCD 
analysis of the corrections to nonleptonic decays of b, is somewhat simplified^ for the ratio 
T n i/T s i, as will be discussed in Sect. 2. This simplification allows to ignore the non-Abelian 
nature of the QCD interaction and to relate the problem to the known solution of a Coulomb 
problem for non-relativistic as well as for relativistic fermions. In Sect. 3 the relativistic case 
is discussed, relevant for the decay b — > cud(s) and in Sect. 4 the decay b — > ccs is considered, 
where both the relativistic and non-relativistic Coulomb dynamics comes into play. In Sect. 5 
the results of this paper are summarized and discussed. 



2 Feynman graphs for T n i/T s i 

The structure of the graphs, which need to be calculated in order to find the QCD corrections 
to T n i/T s i up to 0(a 2 ), becomes transparent if one represents the decay rate as the imaginary 
part of the b quark self-energy. Then the zeroth order nonleptonic decay rate is given by the 
graph of Fig.l, while the 0(a s ) and the 0(a 2 ) corrections are represented by the classes of 
graphs shown respectively in Figures 2 and 30. It is convenient for the purpose of discussion 
to call the quark line b — > c — > b in those graphs as the heavy quark line, and the loop, made 
of ud(s) or cs quarks as the light quark loop. The graphs, where the gluons are attached only 
to the heavy quark line (like those shown in Fig. 2a and Fig.3a), are identical to the same 
graphs for the semileptonic decays and thus they cancel in the ratio T n i/T s i. Therefore the 
corrections to this ratio should contain either gluon exchange within the light quark loop, or 
a gluon exchange between the loop and the heavy quark line, or both. Notice however that 
a color trace is taken over the light quark loop, thus the graphs where an overall color goes 
into the loop do not contribute to the decay rate (examples of such graphs are in Fig. 2c and 
also in Fig. 4). Therefore in the first order in a s only the graph of Fig.2b is left with the gluon 
exchange within the loop. For light quarks this exchange is very well known to result in the 
factor (1 + a s /7r), which enters eq. ([!]), while for the cs loop the charm quark mass enhances 
the correction, so that for realistic quark masses the correction factor is@' 0' (1 + 5 a s /n) 
with 5 « 4.5. 

2 We neglect here the penguin type contributions, originating from the process b — ► ccs, since these are 
known and small § . 
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In the order a 2 the graphs for the corrections to the ratio T n i/Y s i fall into two categories: 
QCD corrections in the order a 2 entirely contained within the loop and a two-gluon exchange 
between the heavy quark line and the loop. Due to the color trace over the loop there is no 
interference in this order of the gluon exchange between the heavy quark line and the loop 
and either the corrections inside the loop or on the heavy quark line. Also the gluon self- 
energy, or non-Abelian splitting of the gluon into two for the gluon exchanged between the 
heavy line and the loop does not contribute due to the color trace. Thus the only exchange 
graphs that are relevant are those with two gluons, where both gluons originate on the heavy 
quark line and end on the quarks within the loop. An example of such graph is shown in 
Fig. 3c. 

The second order QCD corrections contained within the loop are well known for the 
massless quarksE' 0. They replace the first-order factor (1 + a s /7x) by 



with the number of light quark flavors rif = 3 and q 2 being the square of the invariant 
mass of the quarks in the loop. Thus in this case the effect of the second order basically 
reduces to specifying the normalization scale for the coupling a s in the first-order term and 
does not contain any spectacular features beyond that. For the cs loop, where the mass of 
the charmed quark cannot be neglected, the corresponding result for the QCD correction in 
order a 2 is not known. 

Proceeding to the discussion of the graphs with two-gluon exchange between the heavy 
quark line and the light quark loop, we first notice that such graphs have no non-Abelian 
structure and their contribution is exactly the same as it would be in an Abelian theory. The 
only feature of the color structure of QCD that is relevant is the color factor originating from 
the color traces. One can readily see that to correctly reproduce the effect of the two-gluon 
exchange between the heavy line and the loop relative to the zeroth order contribution of 
the diagram of Fig.l, one can replace the QCD interaction by an Abelian theory with the 
coupling a related to a s as a 2 = 2a s /9. (This takes into account the color factor of 3 present 
in the graph of Fig.l). 

A part of the contribution of these graphs is contained in the a 2 term in eq.(p]), that is 
the part, which is due to the exchange of virtual gluons with momenta between rrib and m\y- 
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Here we will be concerned with the opposite part of the spectrum, i.e. with gluons with 
momenta much less than m c . In general, the imaginary part of graphs of the type shown in 
Fig. 3c is contributed by exchange of virtual gluons and by radiation and absorption of real 
transversal gluons. However, for a slow charmed quark its radiation and absorption arises 
only in the order v 2 of the expansion in its velocity^. Thus, the leading terms are determined 
only by the exchange of virtual gluons. Furthermore, it is only the spin-independent (for 
heavy quark) exchange of Coulomb gluons that contributes in the leading order for a slow 
quark. Therefore the problem reduces to the one-particle problem of motion of a quark in 
an Abelian Coulomb-like field of a heavy quark. For the non-relativistic motion of the two 
charmed quarks in the decay b — > ccs this is the standard Quantum Mechanical problem, 
while for the light quarks both in this decay and in b — > cud(s) the problem reduces to the 
known solution of the Dirac equation in a Coulomb field (see e.g. in the textbook |[L4|). In 
the former case the singularity in the velocity of the charmed quarks has the well known form 
n 2 a 2 /v 2 , while in the latter case a logarithmic singularity arises due to the corresponding 
singular behavior of the relativistic Coulomb wave function at the origin in the order a 2 . 



3 The decay b —> cud(s 



We first consider the decay b —>■ cud, where both quarks recoiling against the charmed quark 

are genuinely light (the decay b — > cus is obviously similar in as much as the masses of both 

the d and s quarks can be neglected). As described above, this decay is considered in the 

limit of slow charm, i.e. in the limit A = m& — m c <C m&, m c . In this limit the final charmed 

quark is considered to be at rest, and, as it will be shown, the a 2 term in the correction 

to the decay rate develops a coefficient proportional to ln(m c /A). Certainly, the decay rate 

itself goes to zero in the limit A — > due to the overall A 5 factor, and it is only the ratio of 

the a 2 correction term to the zeroth order decay rate that displays a logarithmic singularity 

in this limit. Also as explained above, in order to find the soft gluon contribution in the 

two-gluon exchange graphs the QCD interaction can be replaced by an Abelian one with 

the coupling a 2 = 2a 2 /9. The interaction of the final light quarks with the charmed quark, 

acting as a Coulomb center, is taken into account by replacing the plane wave functions for 

the light quarks by their wave functions in the Coulomb field. Thus we briefly discuss an 

adaptation of the known textbook results" to the case considered here. Before proceeding 

3 Notice, that the graphs where a gluon is emitted and absorbed by a light quark are not of the type of 
Fig. 3c. In the graphs of this type at least one emission or absorption has to be by a heavy quark. 
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to this discussion, let us note that in the calculation of the Coulomb wave functions of 
the light quarks the charmed quark can be considered as a static center located at r = 0, 
only for the wave functions at the distances r > ro ~ m c~ 1 -, since at shorter distances the 
momenta of the exchanged gluons are not small in comparison with the mass of the charmed 
quark, and its recoil cannot be neglected. Thus we will use the light quark wave functions 
with the logarithmic accuracy down to these distances in the weak interaction Hamiltonian 
proportional to (u'(ro) (1 — 75) d(r )). 

The Dirac-Weyl equation for a massless two-component left-handed d quark spinor in the 
repulsive Coulomb field U(r) —a/r reads as 

{e d -U(r)-i(T- V) d(r) = , (3) 

where Ed is the energy of the d quark. A non-vanishing near the origin solution to this 
equation arises only in the lowest partial wave, where the ansatz for the spinor structure of 
the solution has the form 

d(r) = (f(r) + % g(r) a • n) w (4) 

with w being a constant two-component spinor, and n = r/r. Upon substitution of this 
ansatz in eq.(|^) the resulting equations for the functions f(r) and g(r) exactly coincide with 
the textbook onesG-i in the massless case m = and the angular moment parameter k = — 1. 
Thus the solution can be readily read off the textbook: 

= 25e~- \L n P ' {6^^x^ 7 + ^,27 + 1,-2^ } , 5 

g j T(27 + 1) r R e 1 J 

where 7 = y/1 — a 2 , p = e^, and the phase angle £ is defined as exp(— 2i£) = —7 — ia. The 
wave function of the right handed u antiquark is obtained from this solution by a formal 
replacement a — > —a, g — ► —g. 

In calculating the matrix element of the weak interaction Hamiltonian we are interested 
in these functions at a distance r such that pr 1, but still pr 3> a, since we are seeking 
an expansion in powers of a and if the latter condition were not satisfied such an expansion 
would be impossible. Due to the former condition the degenerate hypergeometric function 
can be replaced by 1. Multiplying the solutions to the wave equations in the expression 
(u*(ro)d(ro)), one finds the ratio of the Coulomb-corrected matrix element to the bare one 
in the following form 



(u^(r )d(r ) 



Coulomb 



(ut(O)d(O)) 



free 



2T( 7 + l + ict 1 



T(2 7 + 1) 



{Ae u e d rly- 1 . (6) 
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One can readily see that the parameter for the expansion of this expression is a 2 rather than 
(a/71") 2 . Under the assumptions made we are however bound to retain only the logarithmic 
part of this expression and only in the order a 2 . With the logarithmic accuracy one can set 
ro ~ wj 1 , £ u ~ £d ~ A thus finding for the correction to the rate of the decay, given by the 
square of the ratio in eq. 



in the form of the factor 
4 



ud 



ni, 



1 + 9 a - hl A 



(7) 



where the relation a 2 = 2a 2 /9 is used. This formula is the final result for the logarithmic 
singularity of the a 2 terms in the decay b — > ud(s) in the slow charm limit of A — > 0. 

Clearly, for the realistic b quark decay the equation (0) can serve only as an indicator 
of possible presence of large terms of order a 2 , rather than as a quantitative estimate of 
these terms, since in reality the charm in b decay is not sufficiently slow for considering 
ln(m c /A) as a large parameter. It is interesting to note however, that the equation (§) can 
be used to find the exact expression for the a 2 correction to a heavy quark decay in an 
entirely artificial kinematical situation, where both the b and c quarks are heavier than the 
W boson: A -C m w -C m b , m c [|. In this limit the W boson propagator provides a cutoff at 
short distances r ~ rn^r, at which the recoil of the heavy quark can still be ignored, thus 
the one-particle wave functions of the form in eq.(|5]) can be used for a complete calculation 
of the matrix element of the weak Hamiltonian. This matrix element is proportional to 



J 



y(r)d(r)) 



exp(-m w r) 3 



d r . 



Calculating this integral with the solutions described by eq.(Q), one readily finds for the 
square of the ratio of the Coulomb-corrected integral to the free one the following expression 



J\£-ui £-d)coulomb 



J(£u, £d)fr 



r(2 7 ) 5 



2T(j + l + ia) 



, m w , 



27-2 



1 + 



A 2 



m 



w , 



(9) 



r(2 7 + i) 

Since the phase space integral with the free wave functions goes in the considered kinematical 
arrangement as e 2 u e\ dsd with e u = A — 6d, one finds the Coulomb correction factor for the 
decay rate as 



T 



r(27)2 3or ( 2 7 + i) 2 



r(4 7 + 2) 



2T( 7 + 1 + ia) 



120 



7 2 r(4 7 + 2) 



|T( 7 + ia) 



r(2 7 + 1) 



/4A 2 



'4A 2, 



27-2 



27-2 



1 2 



(10) 



w , 



4 Actually the condition A <§; mw can be dropped. It is assumed here in order to simplify the resulting 
expression, which anyway serves only for the purpose of illustration. 
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(the identity j 2 + a 2 = 1 is used in the last transition). When expanded to the order a 2 , 
this factor takes the form 

/ 167 7r 2 mw 



30 




21n^ )W S . (11) 



This expression describes both the logarithmic term in the a 2 correction and the non- 
logarithmic one, which is justified to be retained in the assumed artificial limit A -C my/ -C 
77ib, tti c . When combined with the result in eq.@ for the QCD corrections contained within 
the light quark loop, this would completely describe the QCD corrections up to the order a 2 
in this limit. It can be also mentioned that in this limit there of course are no terms with 



4 The decay b —> ccs 

In the decay b — > ccs in the slow charm limit one has two slow quarks in the final state. 
An exchange of gluons between them results in the well-known non-relativistic Coulomb 
corrections, which have as their parameter na s /v, where v is the velocity of either of the 
charmed quarks in their center of mass frame: v = yjl — Am 2 /q 2 with q 2 being the invariant 
mass squared of the cc pair. Thus in the limit of small v the dominant a 2 term is given by 
the double Coulomb exchange between the c and c, and is described by the second term of 
the expansion of the Coulomb factor 

ira/v na n 2 a 2 . . 

/ _, =1 + — + -— + ■■■ (12) 

1 — exp(— ira/v) v 12 v z 

The linear term in a should be discarded, as explained in Sect. 2, while the second term gives 
the correction factor 

1 + n 2 a 2 /(5Av 2 ) (13) 

for the rate (the relation a 2 = 2a 2 /9 is again used here, which takes into account the color 
factors in the relative correction). 

However, one in fact can estimate at least one more term in the expansion of the a 2 
correction in powers of v, and thus get some feeling of what v is sufficiently small for appli- 
cability of the expansion. This is especially relevant, given rather small numerical coefficient 
in the leading correction term in eq.(0). The next term, linear in 1/v arises through the 
interference of the one Coulomb gluon exchange between c and c, described by the term 
0(a) in eq. (|12"D with the less singular than 1/v part of the gluon exchange of either the c 
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or s quarks with the heavy quark line. For the c this corresponds to a hard gluon exchange 
with either the c or the b quark, while for the s quark this corresponds to contribution 
of gluons with generally arbitrary momenta. The contribution of the soft gluons for the s 
quark exchange with the heavy quark line is however not singular at all in v: the logarithmic 
singularity develops only starting from the order a 2 , as described by the wave function in 
eq.(||). Therefore the linear in 1/v term has no additional logarithmic dependence on v. 
The non-logarithmic coefficient can be estimated in the same approximation as in eq.(|l]) i.e. 
using L as a parameter. In this approximation the exchange of the non-Coulomb gluon is 
dominated by the region of momenta between and mw, and one finds that the interfer- 
ence of the hard gluon exchange with a Coulomb exchange within the cc pair adds a negative 
contribution to the correction factor (1131): —3 — L = — § — L. Thus the final result for the 

\| — V v 3 v 

a 2 correction factor to the b — > ccs decay rate, unaccounted for by eq.([L]), in the slow charm 
limit is estimated as 

^ = 1 + (^-£ (i + 0(1)) ) a " • (14) 

The formula in eq. (|T4l ) describes the leading a 2 terms in the part of the spectrum, where 
the cc pair is close to its threshold. In the total decay rate of the real b quark the dominance 
of the singular in v terms is however quite smeared. Indeed, for the average over the spectrum 
of the b — > ccs decay values of v~ 2 and v' 1 one finds at m c /mb ~ 0.3: 

(v~ 2 ) w 5.14 , (v~ l ) « 1.97 . (15) 

Thus the linear in 1/v term in eq. flHD in fact dominates over the quadratic one and makes 
the overall coefficient of a 2 negative and large: F 5s ml — 2.8 a 2 . Certainly, under these 
circumstances one can not rely on the expansion in v for the total rate and may only consider 
eq.(|i~4D as an indication of presence of large contributions in the order a 2 . 

If the coefficients of the a 2 terms may be as big as indicated by the present analysis, these 
terms may well compete with the first-order term in a s . Thus the a 2 effects may essentially 
modify the conclusion based on the first order in a, calculation about the relative 

enhancement of the b — > ccs decay. 

It should be also noted in connection with the b — > ccs decay that the gluon exchange 
contained within the cs loop receives a logarithmic contribution near the threshold in the 
order a 2 from the Coulomb exchange. This contribution also can be readily found from 
eq.(|5]). However in this case one has a competing logarithmic effect due to the so-called 
hybriddll anomalous dimension of the current (s^ 7 m c) both in the leading order and in 
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the next-to-leading ordeiQ. In any event, however, these terms are less singular in the slow 
charm limit than those in eq. (^4|) . 

5 Summary and discussion 

Lacking a complete calculation of the second-order QCD corrections to the b quark decays, 
we have to rely on parametric approximations for these corrections. Unfortunately, in the 
real situation there is no such parameter, which would be acceptably good. The widely used 
expansion of these corrections in L = ln(mw / nib) ~ 2.8 by means of RG does not seem 
to work well, as illustrated by eq. ([!]). In this paper a different parameter is considered: 
the velocity of the charmed quarks in the final state of the decay. The terms which are 
parametrically large in the limit of large 1/v are in general independent of those, which 
dominate the expansion in L, thus they indicate the magnitude of the contributions lost in 
the so-far considered order of the expansion in L. Also in the decay b — > ccs the inverse 
velocity of the charmed quarks is almost as good (or as bad) a parameter as L. As is shown 
the leading in the slow charm limit the a 2 corrections are associated with the interaction 
of the final state quarks with the Coulomb-like gluonic field of the slow quark. Therefore 
the parameter for these corrections turns out to be a 2 in the case of the decay b — > cud(s) 
with relativistic light quarks (eq.(|7])) or (7ra; s ) 2 for the case of the interaction between the 
non-relativistic c and c in the decay b — > ccs (eq.(|H|)). This can be considered as a strong 
indication that the missing in eq. ([TJ) terms of order a 2 can have large coefficients, which may 
be considerably different in the b — > cud(s) decay from those in the b — > ccs decay. In this 
case the missing contributions may well compete in their numerical significance with the 
terms present in eq.(|l|) and may considerably alter the theoretical predictions both for the 
semileptonic branching ratio and for the average charm yield in B decays. Thus the main 
conclusion from the arguments presented in this paper is that it is rather premature to seek 
a contradiction of the theoretical predictions with the data on B s i and n c before a complete 
theoretical analysis of the a 2 corrections becomes available. 

This work is supported, in part, by the DOE grant DE-AC02-83ER40105. 
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Figure 1: The lowest order graph, whose unitary cut describes the rate of the decay b — > cq 2 q\. 
The small filled circles represent the W boson propagators. 
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Figure 2: Three types of graphs, whose unitary cuts describe the first QCD radiative cor- 
rections to the inclusive decay rate b — > cq 2 qi- The dashed lines correspond to gluons. The 
gluon vertices can be anywhere on the be line (a), quark lines in the loop (b), or one vertex 
anywhere on the be line and the other vertex on either line in the loop (c). The graphs of 
the type c in fact do not contribute to the decay rate, as explained in the text. 
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Figure 3: Three types of graphs, whose unitary cuts contribute to 0(a 2 s ) QCD radiative 
corrections to the inclusive decay rate b — > cq 2 Qi- 
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Figure 4: One of the types of graphs of order a 2 , which do not contribute to the decay rate 
of b — > cq 2 qi, because an overall color flows into the colorless quark loop. 
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